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Thermal Structure Analyses for CSM Testbed (COMET) 


This document is the final report for the project entitled "Thermal Structure 
Analyses for CSM Testbed (COMET)," for the period of May 15, 1992 - August 15, 1994. 
The project was focused on the investigation and development of finite element analysis 
capability of the computational Structural Mechanics (CSM) Testbed (COMET) software 
system in the field of thermal structural responses. The stages of this project consisted of 
investigating present capabilities, developing new functions, analysis demonstrations and 
research topics. The Appendices of this report listed the detailed documents of major 
accomplishments and demonstration runstreams for future references. 


1. Summary of Accomplishments 

* Implemented thermal stress analysis capability for solid element ES3 

* Confirmed the thermal stress analysis capability of shell elements ESI, ES5 and ES7 by 
comparing results with classical solution or EAL results. Table I shows the thermal stress 
analysis capability in COMET ESi elements. 


Table l 
ESI 

ES2 

ES3 

ES4 

ES5 

ES6 

ES7 

ES8 

ES10 

Shell 

Shell 

Solid 

Shell 

Shell 

Beam 

Shell 

Shell 

Solid 

Yes 

No 

Yes 

No 

Yes 

Yes 

Yes 

No 

No 


* Confirmed thermal stress analysis of composites in COMET. 

* Demonstrated temperature program input (Appendix I). 

* Completed temperature input from PATRAN to COMET by adding a subroutine 
GCPTEMP in PT2T. All kinds of temperatures which can be generated in PATRAN can 
be transferred to COMET (Appendix II). 

* Developed innovative PATRAN temperature input for temperature varying across the 
thickness (multi-values of temperature at one node.) All kinds of multi-values of 
temperatures of COMET formats can be generated in PATRAN then transferred to COMET 

runstream (Appendix II). 

* Integrated thermal / structural analyses. Using different elements and meshes in thermal 
and structural analyses according to the different thermal and structural loading and 
boundary conditions (Appendix HI). 
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* Developed temperature transformation between different meshes (Appendix HI). 

2. Some Existing Problems 

1 ) The incomplete thermal stress analysis capability for elements ES2, ES4, ES8 and 
ES10. 

2) The incompleteness of PATRAN input. 

PATRAN is a powerful tool for the pre-processing of finite element analysis specially for 
complex structures and applied loads. It has been widely accepted by industry and research 
institutions. If we have a well completed PATRAN interface, it can be expected that more 
users will use COMET. By using PATRAN input user can obtain a COMET runstream, 
this feature greatly helps COMET users. However, the translator PT2T has not been well 
completed, users still need to do some modifications (specially in material properties) on 
the runstream produced from PT2T. 

3) Need of a mini-manual. 

COMET needs a mini-manual to lead new users to quickly start using COMET. 
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Appendix I Programmed Temperatures 


This appendix shows an example of programmed temperature input using GCP 
processor. A rectangular isotropic plate is subjected to a combined mechanical and thermal 
loads as shown in Fig. 1. The results compared well with EAL results and classical 


solutions. 



© 


Boundary conditions: u— v=0 at edges 1 and 2 

Temperature and Me chanical loads; 

T=x+y 

Fx= 0.76923+0.76923x Fy=-0.76923+0.76923x 

Fx=- 1 .53850+0.76923y Fy= 1 .53850+0.76923y 

Fcx=0.7 69231, Fcy= 1 . 538462 

T=x 2 +y 2 

Fx= 0.0 Fy=-0.0027778+3.3333x 2 

Fx=-0.0027778-3.3333y 2 Fy= 0.0 

Fcx=- 1.7542, Fcy=6.5014 

Fig. 1 Rectangular plate subjected combined thermal and mechanical loading 


Case 1 

Edge 4: 
Edge 3: 
NodeC 

Case 2 

Edge 4: 
Edge 3: 
NodeC 
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* COMET runstream 


(programmed temperature load) 
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• rr 


cp /c»m/prc/proclib.gal 
chmod +w proclib.gal 
ima testbed « \andinput 
open/new 1, casel.DBC 
•opan 28, 'proclib.gal' 
•sat plib = 28 
•set echo, off 
*def/i ES_GCP=<true> 

•add recl.clp 

. *set echo on, ma, rod 

•call XI 

•eof 

endinput 



•procedure XI *. 

. Initial Modal Daf. Procedure 

•call ES (function * 'DEFINE ELEMENTS'; es_proc = ESI; es_name = EX47) 


(xqt TAB 

START 231 3 4 5 6 
JLOC 

10 . 0 . 0 . 0 . 1 . 0.11121 
11 2 . 0 . 0 . 2 . 1 . 0 . 


"(xqt GCP 

fabricat ion 
shell 

fabid * 1 
nlayers = 1 
mat id = 1 
thick * 1.0 
intpts * 1 
angle = 0.0 
end 

endfab 

material 

isoal 111 
1.0E7 0.3 0.0 
endmat 


1.0E-6 0.0 0.0 0.0 


(xqt ELD 
<ES_EXPE_CMD> 

NSECT * 1 
1 12 13 2 1 20 10 

(xqt TAB 
CON 1 

ZERO 1 2: 1,11: 12,221,11 
•set echo on, ma, md 
(xqt ESI 

•def/a es_name = EX 4 7 
•do $k *1, 2 

•daf/i es_load_sat = <$k> 

DEFINE TEMPERATURE 

•find 1, DEF.* /seq»dsn 
•ram DEF.* : dsn=<dsn> 

♦find 1, JLOC . STAB . 2 . 5 /seq=dsnjloc 
•ram JLOC . BTAB .2.5 : dsnjloc=<dsnjloc> 

♦do $j»l, 200 200 Elements 

ELEMENT * <$j> 

♦do $i = 12, 15 Node #1 - #4 in DEF.Ex47. 

*def/i ioffl * <$i> 

*gal2mac / Name =nn ode /Iof f =<iof f 1> -- 
/Maxnal 1, <dsn> DATA. <$ j> 

NODE « «$i>-ll> 

*def/i ioff2 = ««nnode>-l>*3» 

*gal2mac /Name=coord(l :2] /Iof f =<iof f 2> /Maxn=2 1, <dsnjloc> -- 
DATA . 1 

*def/f x=<coord(l]> 

♦daf /f y=<coord(2]> 

•if «$k> /eq 1 > /then 

T a «x>*<y» • T=x+y 

•else 

T = <<<x>*<x»-«y>*<y»> 

•endif 

♦enddo 

•enddo 

END_DEF I NE_TEMPERATURE 
•enddo 


[xqt AUS 

SYSVEC : APPL FORC 1 
CASE 1 


•do $i=l, 19 

J = <ll + <ll*<$i>» : <.76923«-<.07o923*<$l>>>, <- . . ■o923'j3* -- 

<.07 6923 *<$ i»> 

•enddo 
•do Si=l , 9 

J=<22i+<Si>>: <-l. 5333 01+<. 2 “6923 *<$!>>>, 
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< . 076923*<$i>» 

•enddo 

j s 231: .769231, 1.538462 

SYSVEC: APPL FORC 2 
CASE 1 
1=2 

•do $i=l, 19 

•def/f xx=<<<$i>*<Si>>*0 . 01> 

J=<ll*ll*<$i»: <0.0027778+<3. 33 33 *<xx>>> 

•enddo 
1 = 1 

•do $i=l, 9 

•def/f yy=«<$i>*<Si»*0 . 01> 

J=<22U<$i»: <0. 0027778 -<3. 3333*<yy>>> 

•enddo 

1=1 2: J=231; -1.7542, 6.5014 

stop 

•call Ii_STATIC ( location=centroids ; nval_meth=l ; stress=<true>; 
print=<true>; load_set=l) 

•call L_STATIC ( location=centroids ; nval_meth=l; stress=<true>; 
print=<true> ; load_set=2) 


•fopen 8 result. all 
•set unit prt=8 
•remark 

•remark RESULTS 

•remark 

•do $k =1, 2 

[xqt VPRT 

lines 70 

•remark 

•remark 

•remark DIS PLACEMENS (U=l, V=2) 

•if «$k> /eq 1 > /then 

•remark (TEMPERATURE DISTRIBUTION T (X, Y) =X+Y AND LOAD CASE 1) 


•else 

•remark (TEMPERATURE DISTRIBUTION T(X, Y) =X*X-Y*Y AND LOAD CASE 
•endif 
•remark 

jointk 2,222,22: 6,226,22: 9,229,22: 11,231,22 
print STAT DISP <$k> 


•remark 

•remark 


(xqt PESR 

select /line*100 
♦remark 


•remark 

•remark STRESSES AT ELEMENT CENTERS 

♦if «$k> /eq 1 > /then 

•remark (TEMPERATURE DISTRIBUTION T(X,Y)=X+Y AND LOAD CASE 


1 ) 


•else 

•remark 

•endif 


(TEMPERATURE DISTRIBUTION T(X, Y) =X*X-Y*Y AND LOAD CASE 


► y rk 

print /id=l: 19:2, 20. 61:79: 2, 80, 121:139:2,140,181:199:2.200 
/load_set=<$k> 


stop 
•enddo 
•fclose 8 
•set unit prt=6 
•end 


2 ) 


2 ) 


ORIGINAL PAGE IS 
Of POOR QUALITY 
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Appendix II PATRAN Temperature Input 


This appendix states the update of PATRAN-COMET temperature transformation 
capability of processor PT2T. As reported before, the processor PT2T has been improved 
to be able to transfer temperature distribution from PATRAN to COMET. Any kinds of 
temperature distribution which can be generated in PATRAN (into neutral files) by using 
"NODE TEMPERATURES" can be transformed into COMET. But PATRAN does not 
generate temperature across the thickness for shell or beam elements, because such 
temperature distribution requires more then one temperature value at one element node and 
PATRAN can only generate one temperature value at one node. 

In the thermal stress analysis of COMET, the temperature distributions across the 
thickness of shell or beam elements are inputted in various formats using DEFINE 
TEMPERATURE. These formats include: (see section HI) 


CONSTANT one value at one node, 

LINEAR two or three values at one node, 

QUADRATIC three values at one node, 

CUBIC four values at one node, 

EXPONENTIAL — three values at one node, 

PWLIN NL+1 values at one node, and 

TABLE NT*2 values at one node 


To transfer such kinds of temperature distributions from PATRAN to COMET, an alternate 
transformation has been developed and the subroutine GCPTEMP in PT2T has been 
further updated to have the capability to transfer all above formats for the temperature 
distribution across the thickness of shell and beam elements form PATRAN to COMET. 
To apply these temperature transformation certain temperature input formats have to be 
followed in PATRAN temperature generation. The basic idea of this new temperature 
transformation is to define different sets of temperatures in PATRAN using designed 
formats. Then to combine them in the translator processor PT2T and obtain the COMET 
runstream for the temperature distribution across the thickness. 

The detailed steps of how to use PT2T to transfer temperature from PATRAN to 
COMET are presented in section H.2. The updated processor PT2T does not affect the 
original temperature distribution transformation. If in a structure temperature does not 
varying across the thickness, user does not need to use these designed formats to input 
temperatures and just generates temperatures using PATRAN formats. The feature of this 
new PATRAN temperature transformation is stated in the part a) of section II.3. The 
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running steps of thermal stress analysis on COMET using PATRAN input is summarized 
in the part b) of section II.3. 

Two examples, a temperature distribution without variation across the thickness is 
given in section EL4 and a temperature distribution with different formats of variation 
across the thickness is given in section H5. 

The updated subroutine GCPTEMP is listed in section n.6. 
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II.l 


COMET Temperature Input Formats 
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User Input for Thermal Stress 


Thermal Stresses 


Table A.l Temperature Input Formats 

Format 

Solid 

Layered Shell 

Beam 

CONST* 

f 

n/ 

f 

V 

r 

V 

LINEAR 

— 



V 

r.Q.G 

V 

PWLIN 

— 


Ti.r, r* t+l 

s/ 

— 


QUADRATIC 

— 


r.,r,,r, 

V 

— 


CUBIC 

— 


rs.Tj.Tj.r, 


— 


EXPON 

— 


«»*.« (T{x) = 


— 


TABLE 

— 


Ti,*!, Tt,x j, Tnt**nt 

v' 

— 



Ker Meaning 

• Default qualifier 

V' Currently implemented 

T Average (ref. surface) temperature 

Tto r Top surface temperature 

7|«t Bottom surface temperature 

N L Number of layers 

NT Number of tabulated points 


A. 2- 2 Generic Constitutive Processor Manual 5/12/1991 
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Therm*! Stresses 


User Input for Therm*] Stress 


For the /LINEAR form.t, lh« .hell Ihrough-thickn.M tempcr.ture di.lribution is 


computed *s follows: 


Tiop + Tut , _ Tt,p ~ Tm 


r(j ) = i i^ + S 

while the beam through-thickness temperature distribution is computed as 

- &T 6T 
r(y.O = ? + + 

For the /P WLIN format, the shell temperatures are interpolated linearly within each 
layer, using the layer interface temperatures, r,(* - l,NL); 

For the /QUADRATIC format, the sheU through-thickness temperature distribu- 
tion is computed as: 

r(») = t, + r, « + r, »* 

For the /CUBIC Iorm»t, the shell through-thickness temperature distribution is 
computed as: 

T{x) ■ T f + T x t + Tit 1 + Is** 


For the /EXPONENTIAL format, the shell through-thickness temperature distri- 
bution is computed via: 

T(z) = • 

Finally, for the /TABLE format, the shell temperatures are tabulated at the through- 

thickness coordinates, z lf z a z/vT, **d piecewise linearly interpolated between these 

points. 

Note that the thickness coordinate z, appearing in the above definitions, is measured 
relative to the shell mid-surface, and h is the shell thickness. 

The structural temperature distribution may be defined with the DEFINE TEM- 
PERATURE command at any time after the element connectivity has been defined 


5 / 12/1991 


Generic Constitutive Processor Manual 


A. 2- 3 
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11.2 

Temperature transformation 
from PATRAN to COMET 
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n Node tem perature variations QnlY 

The word M node temperature variations only" refers: (a) solid elements have temperatures 
defined on nodes only, or (b) sheU or beam elements have temperatures defined in inplane nodes 

only. That means each node has only one temperature value. 

There is no special requirements for these temperature generations, user can just follow 
PATRAN manual to generate temperature distributions. Any kind of temperature distribution 
which can be generated (into neutral file) in PATRAN by using "NODE TEMPERATURES" can 
be transformed into COMET runstream. The new subroutine in PT2T also transfers different 
temperature load sets defined in PATRAN. The only restriction is. 

The number of temperature load set < 10000. 

It is recommended to choose load set number starting with 1, 2, 3, ... sequentially. 

?/) Inplane and across the thickness temperature variations for Shell Of foam 

Temperate, which varies not only in inplane but also across the thickness of shell or beam 
elements, requests to transfer more then one values at one node of the shell or beam elements. 

New subroutine in PT2T transfers those temperatures according to the formats provided in 
COMET (see Attachment I). This transformation requires that user follows designed load set 

numbers to generate temperature load in PATRAN. 

Each designed load set number has five digits. The first digit (lxxxx - 9xxxx) indicates the 
format of the temperature; the second and third digits (xOlxxx - x99xxx) indicate the sequential 
number of input values at one node; the fourth and fifth digits (xxxOl - xxx99) indicate the load set 
number of the temperature distribution. The first, second and third digits are designated and the 
last two digits are user's choice. For example a linear variation of temperature across a shell 
element needs input Ttop and Tbot two values at each nodes. The designated load set numbers for 
PATRAN input are: 

20103 forTtop 

and 20203 for T^ot 

where the first digit 2 denotes the linear variation through the thickness in shell elements; the 
following 01 and 02 are for Ttop and and 03 is load set number in COMET. 

The designed PATRAN load set formats are listed in Table I. 



Table I COMET Temperature Input Formats / Designated PATRAN Load Set Numbers 
COMET Format Element Node Designated load set number 

values for PATRAN input 

( 00 < xx < 100) * 


lOlxx 




CONST 
' = T 


LINEAR 

T(z)= I^3k +z Z^- 


LINEAR 

T(y.z) = T+yf^g 


PWLIN 

T +T T — T. 
T(z) = 1 1 + Z 

hj: Thickness of layer i 


QUADRATIC 

T(z) = Tq+T 1 Z+T 2 Z 2 


♦♦CUBIC 

T(z) = Tq+Tiz+T2Z 2 +T3Z 3 


♦♦EXPON 

T(z) = aeKz<) 


TABLE 


t(z) = Tt+g h+z W T i 
2 


Element 

Node 


values 

SoUd 

Layered Shell 
Beam 

T 

Layered Shell 

Ttop 


Tbot 

Beam 

T 

31 

3y 

31 

dz 

Layered Shell 

Ti 


t 2 


T... 


Tnl+1 

Layered Shell 

To 


Ti 


Layered Shell 


Layered Shell 


Layered Shell 


♦* Not currently implemented in COMET 



801xx 
802xx 
803xx 
804 xx 
8. ..xx 

8(NT*2-l)xx 

8(NT*2)xx (NT<100) 


* xx is COMET load set number 
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In PATRAN input the designated five digit numbers are only for the load set number of 
"NODE TEMPERATURE". They are not the real load set number in COMET. The real load set 
number in COMET are the fourth and fifth digits of the designated five digit number. That means 
the same temperature load set in COMET should have the same fourth and fifth digits in PATRAN 
input For examples, a shell structure has two temperature load sets: Load 1 and Load 5. Load 1 
contains CONST and QUADRATIC temperature distributions and Load 5 has CONST, TABLE 
and LINEAR temperatures. For Load 1 the temperature set 10101 for T, 50101,50201 and 50301 
for TO, T1 and T2 should be generated in PATRAN. And for temperature Load 5 the temperature 
set 10105 for T; 80105, 80205, 80305, 80405, ... for Tl, Zl, T2, Z2,... and 20105, 20205 for 
Ttop,Tbot should be generated in PATRAN. 

Basically this transformation treats node temperature multi-values as different designed 
temperature sets in PATRAN, then combine them in PT2T to become COMET formats. 
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II.3 


Running thermal stress analysis in COMET 
using PATRAN input 
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3) Type 'pt2ttemp' to execute translator PT2T. A runstream file (e.g. PT2T.PRC) will be 
created. 

4) Modify the COMET runsteam in PT2T.PRC or prepare user’s own runstream file which calls 

the transformed subroutines in PT2T.PRC. 

5) Create an executable file xxx.com, and execute an analysis. 
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II. 4 

An example with general temperature distributions 
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b=1 0 



CONST 

T(Z) = T 

LINEAR 

T(z) = T ‘° p ^ Tb ° 1, + z Tl ° P h Tb0 - 

PWLIN 

T(z) = + Z T[+ l T± 

hi: Thickness of layer i 
hi=h3=0.03, h2=0.04 

QUADRATIC 

H 

/—s 

N 

V— ' 

II 

I 

N 

TABLE 

T(Z) 2 
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1) TEMP.ses PATRAN input 
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GO 

1 

1 

2 

PA, 1,, 10/10 
GFEG,P1, ,11/6 
CFEG , P 1 , QUAD 

DFEG,P1,DISP, 0/0/0, ,EDl/4 

PF, PI, QU, 1/0.1 

PMAT, 1 , ISO, 3 • 0E7 , , 0. 3, , 6 . OE-6 

FIELD, 2, RECT/S, 1/100/1/100 
DFEG, PI, TEMP /N,S2, 10101, H1T18 
DFEG, PI, TEMP/N, 200, 20101, N19T30 
* DFEG, PI, TEMP/N, -200, 20201, N19T30 
DFEG, PI, TEMP /N, 410, 40101, N31T42 
DFEG, PI, TEMP/N, 420, 40201, N31T42 
DFEG, PI, TEMP/N, 430, 40301, N31T42 
DFEG, PI, TEMP/N, 440, 40401, N31T42 
DFEG, PI, TEMP/N, 510, 50101, N43T54 
DFEG, PI, TEMP/N, 520, 50201, N43T54 
DFEG, PI, TEMP/N, 530, 50301, N43T54 
DFEG, PI, TEMP/N, 810, 80101, N55T66 
DFEG, PI, TEMP/N, -0.05, 80201, N55T66 
DFEG, PI, TEMP/N, 830, 80301, N55T66 

DFEG, PI, TEMP/N, -0.02, 80401, N55T66 

DFEG, PI, TEMP/N, 850, 80501, N55T66 
DFEG, PI, TEMP/N, 0.01, 80601, N55T66 

DFEG, PI, TEMP/N, 870, 80701, N55T66 
DFEG, PI, TEMP/N, 0.05, 80801, N55T66 
DFEG, PI, TEMP/N, 5 12,5 0102,N43T54 
DFEG, PI, TEMP/N, 52 2, 502 02, N43T54 
DFEG, PI, TEMP/N, 532, 50302, N43T54 
DFEG, PI, TEMP/N, 812, 80102, N55T66 

DFEG, PI, TEMP/N, -0.05, 80202, N55T66 

DFEG, PI, TEMP/N, 832 , 80302 ,N55T66 
DFEG, PI, TEMP/N, -0.02, 80402, N55T66 
DFEG, PI, TEMP/N, 852, 80502, N55T66 
DFEG, PI, TEMP/N, 0.01, 80602, N55T66 
DFEG, PI, TEMP/N, 872, 80702, N55T66 
DFEG, PI, TEMP/N, 0.05, 80802, N55T66 
END 


/ 
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2) TEMP.PRC PT2T output 
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♦PROCEDURE PT2T — START 
START 66 
♦END 

♦PROCEDURE PT2T_TITLE 
TITLE' TEKP1 
♦END 

♦PROCEDURE PT2T_JLOC 


JOINT LOCATIONS 

i 0*0 

0.0 

0.0 

2 

0.0 

2.000000 

0.0 

3 

0.0 

4.000000 

0.0 

4 

o.o 

6.000000 

0.0 

5 

0.0 

8.000000 

0.0 

r':? 6 

0.0 

10.00000 

0.0 

7 

1.000000 

0.0 

0.0 

8 

1.000000 

2.000000 

0.0 

9 

. 1.000000 

4.000000 

0.0 

10 

1.000000 

6.000000 

0.0 

11 

1.000000 

8.000000 

0.0 

12 

1.000000 

10.00000 

0.0 

13 

2.000000 

0.0 

0.0 

f - 14 

2.000000 

2.000000 

0.0 

15 

2.000000 

4.000000 

0.0 

16 

2.000000 

6.000000 

0.0 

17 

2.000000 

8.000000 

0.0 

18 

2.000000 

10.00000 

0 . 0 

19 

3.000000 

0.0 

0.0 

20 

3.000000 

2.000000 

0.0 

21 

3.000000 

4.000000 

0.0 

22 

3.000000 

6.000000 

0 . 0 

23 

3.000000 

8.000000 

0.0 

24 

3.000000 

10.00000 

0.0 

25 

4.000000 

0.0 

0.0 

26 

4.000000 

2.000000 

0.0 

1. ■ 27 

4.000000 

4.000000 

0.0 

r.v 28 

4.000000 

6.000000 

0.0 

29 

4.000000 

8.000000 

0.0 

30 

4.000000 

10.00000 

0.0 

31 

5.000000 

0.0 

0.0 

| 32 

5.000000 

2.000000 

0.0 

33 

5.000000 

4.000000 

0.0 

34 

5.000000 

6.000000 

0.0 

3S 

5.000000 

8.000000 

0.0 

36 

5.000000 

10.00000 

0.0 

37 

6.000000 

0.0 

0.0 

38 

6.000000 

2.000000 

0.0 

39 

6.000000 

4.000000 

0.0 

40 

6.000000 

6.000000 

0.0 

41 

6.000000 

8.000000 

0.0 

42 

6.000000 

10.00000 

0.0 

I 43 

7.000000 

0.0 

0.0 

44 

7.000000 

2.000000 

0.0 

45 

7.000000 

4.000000 

0.0 

46 

7.000000 

6.000000 

0.0 

47 

7.000000 

8.000000 

0.0 

1= . 48 

7.000000 

10.00000 

0.0 

49 

8.000000 

0.0 

0.0 

50 

8.000000 

2.000000 

0.0 

51 

8.000000 

4.000000 

0.0 

i-r c . 52 

8.000000 

6.000000 

0.0 

53 

8.000000 

8.000000 

0.0 

54 

8.000000 

10.00000 

0.0 

55 

9.000000 

0.0 

0.0 

' 56 

9.000000 

2.000000 

0.0 

F .‘ 5 7 

9.000000 

4.000000 

0.0 

58 

9.000000 

6.000000 

0 . 0 
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59 9*000000 8*000000 

(0 9.000000 10.00000 

61 10.00000 0.0 

62 10.00000 2.000000 

63 10.00000 4.000000 

64 10.00000 6.000000 

65 10.00000 8.000000 

66 10.00000 10.00000 

♦END 

♦PROCEDURE PT2T_MAT 
MATERIAL CONSTANTS 

l 3.0000+7 

♦END 

♦PROCEDURE PT2T_TEMP 

[XQT <ES PROO 

♦DEF/A~ ES_NAME - <ES_NAME> 

♦DEE /I ES~LOAD_SET - 1 

DEFINE TEMPERATURE 
ELEMENT - 1 

NODE - 1 

T - 0 . OOOOOOOOE+OO 

NODE - 2 

T • 0. 10000000E+03 

NODE - 3 

T - 0. 30000000E+03 

NODE - 4 

T - 0.20000000E+03 

ELEMENT - 2 

NODE - 1 

T « 0. 10000000E+03 

NODE - 2 

T - 0.20000000E+03 

NODE - 3 

7 - 0. 40000000E+03 

NODE - 4 

T - 0. 30000000E+03 

ELEMENT • 3 

NODE - 1 

T - 0. 20000000E+03 

NODE - 2 

T/LINEAR « 0. 20000000E+03 

NODE - 3 

T/LINEAR - 0. 20000000E+03 

NODE * 4 

T - 0.40000000E+03 

ELEMENT - 4 

NODE - 1 

T/LINEAR ■ 0.20000000E+03 

NODE - 2 

T/LINEAR - 0.20000000E+03 

NODE « 3 

T/LINEAR - 0 . 20000000E+03 

NODE « 4 

T/LINEAR ■ 0 . 20000000E+03 

ELEMENT “ 5 

NODE - 1 

T/LINEAR - 0 . 20000000E+03 

NODE - 2 

T/PWLIN « 0. 41000000E+03 

0. 42000000E+03 
0 . 43000000E+03 
0 . 44000000E+03 

NODE - 3 

T/PWLIN * 0. 41000000E+03 

0. 42000000E+03 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


0.300000 6.0000-6 


-0.20000000E+03 

-0.20000000E+03 


-0.20000000E+03 
-0 . 20000000E+03 
-0 . 20000000E+03 
-0 . 20000000E+03 


-0 . 20000000E+03 

++ 

++ 

++ 

++ 


++ 

++ 
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T/TABLE 


0*810000002+03 -0. 50000001E-01 ++ 
0 . 83000000E+03 - 0 . 20000000E-01 ++ 
0*850000002+03 0.999999982-02 ++ 
0.870000002+03 0.50000001B-01 ++ 


: 

NODE - 4 

T/ QUADRATIC ' 
ELEMENT - 
NODE - 1 

T/TABLE - 


« 0 . 51000000E+03 
10 

0.810000002+03 -0. 
0.830000002+03 -0. 
0. 85000000E+03 0. 
0.870000002+03 0. 


0.520000002+03 


500000012-01 ++ 
20000000E-01 ++ 
99999998E-02 ++ 
50000001E-01 ++ 


NODE - 
T/TABLE * 


2 

0.81000000E+03 -0. 50000001E-01 ++ 
0.830000002+03 -0. 20000000E-01 ++ 
0.850000002+03 0. 99999998E-02 ++ 
0.870000002+03 0.50000001E-01 ++ 


NODE - 
T/TABLE - 


0.810000002+03 - 0 . 50000001E-01 ++ 
0.830000002+03 -0 . 20000000E-01 ++ 
0.850000002+03 0. 99999998E-02 ++ 
0.870000002+03 0 . 50000001E-01 ++ 


0 . 53000000E+03 


NODE - 
T/TABLE - 


0.810000002+03 -0. 50000001E-01 ++ 
0.830000002+03 - 0 . 20000000E-01 ++ 
0.850000002+03 0. 99999998E-02 ++ 
0.870000002+03 0. 50000001E-01 ++ 


T 

T 

T 

T 


ELEMENT - H 

NODE - 1 

• 0 . 20000000E+03 

NODE - 2 

- 0 . 30000000E+03 

NODE - 3 

. 0.500000002+03 

NODE - 4 

a 0. 40000000E+03 
ELEMENT - 12 


NODE - 1 

T a 0. 30000000E+03 
NODE - 2 

T - 0 . 40000000E+03 

NODE - 3 

T - 0 . 60000000E+03 

NODE *• 4 

T - 0 . 50000000E+03 

ELEMENT “ 13 

NODE ■ 1 

T - 0. 40000000E+03 

NODE - 2 

T/LINEAR * 0.20000000E+03 

NODE - 3 

T/LINEAR “ 0 . 20000000E+03 

NODE - 4 

T - 0 . 60000000E+03 

ELEMENT - 14 

NODE - 1 

T/LINEAR ■ 0 . 20000000E+03 

NODE - 2 

T/LINEAR * 0 . 20000000 E+ 03 

NODE - 3 

T/LINEAR “■ 0.200000002+03 


-0.200000002+03 

-0.200000002+03 


-0 . 20000000E+03 
-0.200000002+03 
-0 .200000002+03 
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1 


m 0 • 51200000E+03 0.52200000E+03 0. 53200000E+03 

^0. 81200000E+03 -0.50000001E-01 +♦ 

0.83200000E+03 -0 . 20000000E-01 ++ 

0 . 85200000E+03 0. 99999998E-02 +♦ 

0.87200000E+03 0 . 50000001E-01 ++ 


NODE - 
T/TABLE « 


0 . 81200000E+03 -0 . 50000001E-01 ♦♦ 
0 . 83200000E+03 - 0 . 20000000E-01 +♦ 
0 . 85200000E+03 0 . 99999998E-02 ++ 
0. 87200000E+03 0 . 50000001E-01 ++ 


NODE - 4 

T/QUADRATIC - 
ELEMENT - 
NODE - 1 

T/TABLE - 0 


0.51200000E+03 

50 


0.52200000E+03 0. 53200000E+03 


0. 81200000E+03 -0. 50000001E-01 ++ 
0. 83200000E+03 -0 . 20000000E-01 ♦+ 
0 . 85200000E+03 0. 99999998E-02 ++ 
0 . 87200000E+03 0. 50000001E-01 ++ 


NODE - 
T/TABLE 


0 . 81200000E+03 -0 . 50000001E-01 ++ 
0 . 83200000E+03 -0 . 20000000E-01 ++ 
0 . 85200000E+03 0 . 99 999998E-02 ++ 
0.87200000E+03 0 . 50000001E-01 ++ 


NODE - 
T/TABLE 


0. 81200000E+03 - 0 . 50000001E-01 ++ 
0 . 83200000E+03 -0 . 20000000E-01 ++ 
0 . 85200000E+03 0. 99999998E-02 ++ 
0 . 87200000E+03 0 ♦ 50000001E-01 ♦+ 


NODE - 
T/TABLE 


0 . 81200000E+03 -0 . 50000001E-01 ++ 
0. 83200000E+03 -0. 20000000E-01 ++ 
0 . 85200000E+03 0. 99999998E-02 ++ 
0. 87200000E+03 0. 50000001E-01 ++ 


END_DEFINE_TEMPERATURE 
♦END ~ 

♦PROCEDURE PT2T_CONN 

<ES EXPE CMD> ~~ • QUAD4 


ELEMENT TYPE 


NMAT- 


1 



NSECT- 


1 



1 

7 

8 

2 


7 

13 

14 

8 


13 

19 

20 

14 


19 

25 

26 

20 


25 

31 

32 

26 


31 

37 

38 

32 


37 

43 

44 

38 


43 

49 

50 

44 


49 

55 

56 

50 


55 

61 

62 

56 


2 

8 

9 

3 


8 

14 

15 

9 


14 

20 

21 

15 


20 

26 

27 

21 


26 

32 

33 

27 


32 

38 

39 

33 


38 

44 

45 

39 


44 

50 

51 

45 


50 

56 

57 

51 

• 


PATRAN 

PATRAN 

PATRAN 

PATRAN 

PATRAN 

PATRAN 

PATRAN 

PATRAN 

PATRAN 

PATRAN 

PATRAN 

PATRAN 

PATRAN 

PATRAN 

PATRAN 

PATRAN 

PATRAN 

PATRAN 

PATRAN 
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I 


56 

62 

63 

57 


3 

9 

10 

4 


9 

15 

16 

10 


15 

21 

22 

16 


21 

27 

28 

22 


27 

33 

34 

28 


33 

39 

40 

34 


39 

45 

46 

40 


45 

51 

52 

46 


51 

57 

58 

52 


57 

63 

64 

58 


4 

10 

11 

5 


10 

16 

17 

11 


16 

22 

23 

17 


22 

28 

29 

23 


28 

34 

35 

29 


34 

40 

41 

35 


40 

46 

47 

41 


46 

52 

53 

47 


52 

58 

59 

53 


58 

64 

65 

59 


5 

11 

12 

6 


11 

17 

18 

12 


17 

23 

24 

18 


23 

29 

30 

24 


29 

35 

36 

30 


35 

41 

42 

36 


41 

47 

48 

42 


47 

53 

54 

48 


53 

59 

60 

54 


59 

65 

66 

60 



*END 

* PROCEDURE PT2T_BC 
CONSTRAINT DEFINITION 1 
ZERO 123 
1 
2 

3 

4 

5 

6 
7 

13 

19 

25 

31 

37 

43 

49 

55 

61 

♦END 

♦PROCEDURE PT2T_MODEL 

♦CALL ES ( FUNCTION * 'DEFINE ELEMENTS' ; — 

ES_NAME - <ES_NAME> } — 

ES_PROC * <ES_PROC> ? — 

ES__PARS « <ES~PARS> ) 

[XQT TAB 

♦CALL PT2T_START 
♦CALL PT2T TITLE 
♦CALL PT2T~JLOC 
♦CALL PT2T_BC 
STOP 

♦CALL MAT DAT . USER DEFINED SECTION PROPERTIES. 
[XQT LAU 
[XQT ELD 


PATRAN ELEMENT #20 
PATRAN ELEMENT #21 
PATRAN ELEMENT #22 
PATRAN ELEMENT #23 
PATRAN ELEMENT #24 
PATRAN ELEMENT #25 
PATRAN ELEMENT #26 
PATRAN ELEMENT #27 
PATRAN ELEMENT #28 
PATRAN ELEMENT #29 
PATRAN ELEMENT #30 
PATRAN ELEMENT #31 
PATRAN ELEMENT #32 
PATRAN ELEMENT #33 
PATRAN ELEMENT #34 
PATRAN ELEMENT #35 
PATRAN ELEMENT #36 
PATRAN ELEMENT #37 
PATRAN ELEMENT #38 
PATRAN ELEMENT #39 
PATRAN ELEMENT #40 
PATRAN ELEMENT #41 
PATRAN ELEMENT #42 
PATRAN ELEMENT #43 
PATRAN ELEMENT #44 
PATRAN ELEMENT #45 
PATRAN ELEMENT #46 
PATRAN ELEMENT #47 
PATRAN ELEMENT #48 
PATRAN ELEMENT #49 
PATRAN ELEMENT #50 
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'DEFINE FREEDOMS' ) 


♦CALL PT2T_CONN 
•CALL ES (“FUNCTION - 
(XQT AUS 
•CALL PT2TJTEMP 
•END 
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cp /csm/prc/proclib.gal . 
chmod +w proclib.gal 
time testbed « \endinput 
♦open/new 1, temp.DBC 
♦open 28, 'proclib.gal' 
♦set plib ■ 28 
♦set echo, off 
♦def/i ES_GCP-<true> 

♦add PT2T.PRC 
♦add temp.clp 

. *set echo on, ma, md 

♦call TEMP 

♦eof 

endinput 


32 


♦procedure TEMP 
♦def/a es_jproc - 
♦def/a es_name - 
♦call ES (function 
es name* 


ESI 

EX47 

- 'DEFINE ELEMENTS 
<es name>) 


[xqt TAB 
START 66 
. START 121 
. *call PT2T_START 
♦call PT2TJTITLE 
♦call PT2T__JLOC 
♦call PT2T_BC 


es_proc - 


stop 

. ♦call ES ( function - 
. [xqt ESI 

. *def/a es_name * EX47 
. ♦call MATDAT 
. [ xqt LAU 
[xqt GCP 

fabrication 
shell 
f abid-1 


'DEFINE FREEDOMS ' ) 


nlayers-3 

matid-3*l 

thick-0.03,0.04,0.03 
intpts-3®2 
angle-3® 0. 0 
end 

endfab 
material 
isoel 111 

3.0E7 0.3 0.0 6.0E-6 0.0 0.0 0.0 


endroat 


[xqt ELD 
<ES__EXPE CMD> 

♦call PT2T_CONN 

♦set echo on, ma, md 


♦call PT2T__TEMP 
. [ xqt AUS 
. ♦call PT2T_AD 
. *call PT2T_AF 

♦call L STATIC ( stress-<true>? 


print=<f alse> ) 


3 RECNAME - N0DES_S1 


♦call STRESS 
[xqt NVAL 

reset DIM — 2 METHOD 
DSNAME - STRS.EX47 .♦.* 

T xot PNSR 

1 print /id-5:U5:ll, 56:66:2 /METHOD 
print /id- 3 s 63:6,31:36 /METHOD - 3 


stop 

♦end 


<esjproc> 
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II. 5 

An example with inplane temperature only 


34 



Inplane Tempertature Variation Only 



a = 1.0 
b = 3.0 
T = 100 


1/4 Model 



200r 
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. 300046875E+03 



T • 0. 320051605E+03 
ELEMENT - 5 

NODE • 1 

X - 0.300054260E+03 

X°- E 0 . 300051758E+03 

?°- E 0 . 32 00 5 69 4 6E+ 03 

T°- E 0.320059692E+03 

element • 8 

NODE • 1 

X - 0. 300051758E+03 

NODE • 3 

X - 0.300041534E+03 

NODE - 3 

X - 0.32004S715E+03 

NODE * 4 

X - 0.320056946E+03 

ELEMENT - 1 

NODE - 1 

X - 0.300041534E+03 

NODE • 2 

X - 0 . 300027405E+03 

NODE - 3 

X - 0. 3200301 82E+ 03 

NODE - 4 

X - 0 . 32004S715E+03 
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ELEMENT - 
NODE - 1 , 

X - 0.320051636E+03 
NODE -2 
X - 0 . 320033600E+03 
NODE - 3 

X - 0. 340036682E+03 
NODE * 4 ^ 

X - 0.340056305E+03 

ELEMENT * 39 

NODE - 1 

X - 0. 320033600E+03 

NODE - 2 

X - 0.320011841E+03 

NODE - 3 

X - 0 . 340012939E+03 

NODE - 4 

X - 0. 340036682E+03 

ELEMENT - 40 

NODE - 1 

X - 0. 32 0011 84 1E+ 03 

NODE - 2 

X - 0 . 320000000E+03 

NODE - 3 

T - 0.340000000E+03 

NODE - 4 

X - 0.340012939E+03 


ELEMENT • 51 

NODE - 1 

X - 0- 340000031E+03 

NODE - 2 

X - 0.340004333E+03 

NODE • 3 

T - 0.360004730E+03 

NODE - 4 

T » 0 . 360000031E+03 


ELEMENT ■ "7 2 
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node- • 1 

T - 0. 3600047306+03 
NODE - 2 

X - 0.3600176706+03 

X°- E 0. 3800190126+03 
NODE • 4 

X - 0 . 3800050356+03 
ELEMENT - 7 3 

NODE - 1 

X - 0.360017670E+03 
NODE - 2 ^ 

X - 0.360035706E+03 
NODE “ 3 

X - 0.380038422E+03 
NODE - 4 

X - 0.3800190126+03 


ELEMENT - 94 

NODE - 1 „ 

X - 0.380038422E+03 
NODE - 2 

X - 0.3800581976+03 
NODE “ 3 

X - 0.4000 62 317E+03 
NODE - 4 

X - 0.400041138E+03 


ELEMENT * 295 

NODE - 1 

X - 0 . 580099 60 9E+ 03 

NODE - 2 

X - 0.580124084E+03 

NODE - 3 

X - 0 . 600129333E+03 

NODE - 4 

X - 0.600103821E+03 


ELEMENT - 296 

NODE - 1 

X - 0.580124084E+03 
NODE « 2 

X - 0.580130127E+03 
NODE - 3 _ 

T - 0 . 600135559E+03 
NODE - 4 

T - 0 . 600129333E+03 


ELEMENT - 397 

NODE - 1 

T - 0.680157288E+03 

NODE - 2 

T - 0.680135925E+03 

NODE - 3 

T - 0.7001 4 0686E+03 

NODE - 4 

T - 0.700162781E+03 

ELEMENT - 398 

NODE - 1 

T - 0.680135925E+03 

NODE - 2 

T - 0 . 680088501E+03 

NODE - 3 

T *■ 0.7000915S3E+03 

NODE - 4 

T - 0 . 70014C686E+03 

ELEMENT - 359 
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PATRAS ELEMENT # 111 

116 117 138 137 ; PATRAS ELEMENT # 112 

117 118 139 138 . PATRAS ELEMENT # 113 

118 119 140 139 . PATRAS ELEMENT # 114 

119 120 141 140 • PATRAS ELEMENT #. 115 

120 121 142 141 




.PROCEDURE P T 2 T_BC 
CONSTRAINT DEFINITION 


ZERO 

1 

ZERO 

21 

ZERO 

22 

ZERO 

42 

ZERO 

43 

ZERO 


€3 

ZERO 

64 

ZERO 

84 

ZERO 

85 

ZERO 

105 
ZERO 

106 
ZERO 

126 

ZERO 

127 

ZERO 

147 

ZERO 


2 

1 

2 

1 

2 

1 

2 

1 

2 


148 

ZERO 1 

168 

ZERO 2 

169 

ZERO 1 

189 

ZERO 2 

190 

ZERO 1 

210 

ZERO 2 

211 

ZERO 1 


231 

ZERO 2 

232 

ZERO 1 

252 

ZERO 2 

253 

ZERO 1 

273 

ZERO 2 

274 

ZERO 1 

294 

ZERO 2 

295 

ZERO 1 

315 

ZERO 2 

316 

ZERO 1 

336 

ZERO 2 

337 

ZERO 1 

357 

ZERO 2 

358 

ZERO 1 

378 

ZERO 2 


1 
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1 



/ 


379 ' 
ZERO 

399 
ZERO 

400 
ZERO 

420 
ZERO 

421 
ZERO 
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2 

1 

2 

1 


•END 

* PROCEDURE PT2T__ MODEL 
•CALL ES ( FUNCTION - 
ES NAME “ 
es*”proc * 
es~pars - 


• DEFINE ELEMENTS' 

<es_name> 

<e$_proc> 

<ES*”PARS> 



[XQT TAB 

•CALL PT2T_START 
•CALL PT2T~TITLE 
•CALL PT2T_JL0C 
•CALL PT2T_BC 

!SL MATDAT . USER DEFINED SECTION PROPERTIES. 

[XQT LAU 
[XQT ELD 

•SlL ES^FUNCTION - 'DEFINE FREEDOMS' ) 

[XQT AUS 

•CALL PT2T_TEMP 

•END 
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11.6 

Subroutine GCPTEMP 
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2 


SUBROUTINE GCPTEMP ( TEHPN , NOID, LDATA, NEHS , NGR, NEL, NISID/ IT) 
COMHOM /FXlS/ NINP/NOUT/NPRT/NEUF/NCOH/NVES/NRPT/NLOD/NE43/NUKD/ 

+ NLD1 / NLD2 , NLD3 

COMMON/HEDATA/ ITC, NID, ISID 
COMMON /STEVE/ LFKEY, JBC, JAD, JAF , JAT 
COMMON /CSTEVE/ KFKEY 
CHARACTER KFKEY* 40 

HISIO(NN), STEHPtNNO), LD«A(3..., 
DlJffiNS ION TN02^ 20^32) ^STEMP2 ( 2 0 , NNO ) , J IT ( NN , 3 ) , HAXNO (NNO) 
CAS N CARDOTr ^ JpROCEDURE 1 ™/^^ ( 1 : LFKEY ) // ' _TEMP ' ) 


WRITE (NCOM/ '(A)') ' [XQT <ES_PROO 
WRITE (NCOM, '(A)') ' *def/a ES_NAME 
IF (NN.LT.IT) WRITE(NCOM,'(A,I20)') 

^IF(NNO.LT.NGR) WRITE(NCOM/ ' (A,I20) ' ) 


- <ES NAME> ' 

'NN IS~TOO SMALL, SHOULD > ' 
'NNO IS TOO SMALL, SHOULD >' 


$ ,NGR 
IMN-0 

DO 2 I-1,IT 

GJNN-NISID(I)/100.0 

JNN-INT(GJNN) 

JN3-NISID(I)-JNN*100 

GJN1-JNN/ 100.0 

JNl-INT(GJNl) 

JN2-JNN-JN 1*100 


JIT(I, 1)-JN1 
JIT(I,2)-JN2 
JIT(I,3)-JN3 
IMN-IMN+JN1 
CONTINUE 

IF(IMN.EQ.O) GOTO 1000 


15 

16 

17 


110 


120 


140 


DO 16 1-1, IT 
IF (NISID(I).GT.O) GOTO 17 


GOTO 1100 
JJ* JIT{ I, 3 ) 

WRITE (NCOM, ' (A, 15) ' ) 
WRITB(NCOM, ' (A) ' ) ' 
DO 110 1-1, NGR 
MAXNO(I)-0 
J1NOID(I)-0.0 
DO 110 J-l, 20 


*def/i es_load_set 
DEFINE TEMPERATURE' 


STEMP2 ( J, I ) - 0.0 
DO 120 1-1, IT 
JJ1-JIT(I,1) 

JJ2-JIT(I,2) 

IF(JJ.NE.JIT(I,3) ) GOTO 120 
STEMP2 ( JJ2 ,NOID( I ) ) - TEMPN(I) 
JlNOID(NOID(I) )-JJl 
MAXNO(NOID(I) )-MAXNO(NOID(I) )+l 
NISID(I)«-NISID(I) 


CONTINUE 

DO 130 NE-1, NEL 

NJ1T-0 

DO 140 1-1, NENS(NE) 

J1TN02 ( I )*JlNOID(LDATA( 6+1 ,NE) ) 
N J 1T-N J 1T+ J 1TN02 ( I ) 


CONTINUE 

IF (NJ1T.EQ. 0 ) GOTO 130 
WRITE ( NCOM, '{A, 18)') ' 
DO 145 1*1, NENS(NE) 


ELEMENT - ' , NE 


J*1 


JJ 
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155 

150 

210 

220 

230 

240 

241 

250 

260 

270 

2S0 

281 

290 

145 

130 


IF(J1TN02( I) .EQ.O) GOTO 145 
WRITE(NCOM, ' (A, 13) * ) ' MODE- , I 
TN02 ( J, I ) - STEMP2 ( J ,LDATA( 6+1 ,NE) ) 

IF(J.GE.MAXNO(LDATA(6+I,NE) ) ) GOTO 160 

J-J+l 
GOTO 155 

PRINT 2 *? ' X^ALID^ TEKPRATURE 0 ^! 1 ^ CARD TYM^OOHD ON NEUTRAL FILE' 
GOTO 290 

WRITE(NCOM, ' (A # E16 . 8 ) ' ) ' T * ',TN02(1,I) 

WRITE (NCOM, ' (A,2(El6.8,2x) ) ' ) ' T/LINEAR - ' , (TN02 ( K, I ) /K-l, 2 ) 
WRITE (NCOM, ' (A,3(El6.8,2x) )' ) ' T/LINEAR - ' , (TN02 ( K, I) , K-l , 3 ) 

WRITE(NCOM r ' (A,E16 . 8 , A) ' ) ' T/PWLIN - ' ,TN02( 1, 1) , ' ++' 

DO 241 K-2,MAXNO(LDATA(6+I,NE)) # TW02 ( K I\ ' ++' 

WRITE (NCOM, '(A,E16.8,A)') ,TN02(K,I), 

WRITE (NCOM, * ) ' s' 

WRITE (NCOM, ' (A,3(E16.8,2x) ) ' ) ' T /QUADRATIC - ' , (TN02 (K, I ) , K-l, 3 ) 

WRITE(NCOM, ' (A, 4E16. 8) ' ) ' T/CUBIC - ' , (TN02(K, I ) , K-l , 4 ) 

WRITE(NCOH, ' (A, 3E16 . 8 ) ' ) ' T/EXPON - ' , (TN02<K, I) ,K«1,3) 

WRITE(NCOM, ' (A,2E16.8,A) ' ) ' T / TABLE - ' , (TN02 { K, I ) , K-l , 2 ) , ' ++' 
DO 281 K-3, MAXNO( LDATA ( 6+1, NE) ),2 

WRITE (NOOK, ' (12x,2E16.8,A) ' ) TN02(K,I) ,TN02(K+1,I) , ' ++ 

WRITE(NCOM,*) ' *' 

CONTINUE 

CONTINUE 

C WRITE(NCOM, ' (A) ' ) ' END_DEFINEJTEMPERATURE' 

GOTO 15 


NGR 


JJ 


1100 WRITE(NCOM, ' (A) ' ) ' *END' 

GOTO 1200 
1000 CONTINUE 

5 DO 6 I-l, IT 

6 IF (NISID(I) .GT.0) GOTO 7 
GOTO 100 

7 JJ- NISID(I) 

C WRITE (NCOH, ' (A, 15) ' ) 'NGR- 

WRITE (NCOM, ' (A, 15 ) ' ) ' ‘def/i es load set 
WRITE(NCOH, ' (A) ' ) ' DEFINE TEMPERATURE 

DO 10 I-l, NGR 
10 STEMP(I) - 0.0 
DO 20 I-l, IT 

IF( JJ.NE.NISID(I) ) GOTO 20 
STEMP ( NOID ( I ) ) - TEMPN(I) 

NISID(I)— -NISID(I) 

20 CONTINUE 

DO 30 NE-1, NEL 
TTNO-O . 0 

DO 40 I-l, NENS(NE) 

t* WttTTFfNCOM 9 f A, 13 ) 9 ) # LDATA* 1 ^ LDATA( 6+1 # NE ) 

c TOITE(NCOm''(A,E16.9)' ) 'STEMP (LDATA)-' , STEMP (LDATA ( 6+1 ,NE) ) 

TNO(I) - STEMP ( LDATA ( 6+1 ,NE) ) 

40 TTNO-TTNO+ ABS ( TNO ( I ) ) 

IF ( TTNO . EQ .0.0) GOTO 30 

WRITE (NCOM, ' ( A, 18 ) ' ) ' ELEMENT NE 

DO 45 I-l, NENS(NE) 

IF (TNO( I ) . EQ. 0.0) GOTO 45 
WRITE(NCOM, ' (A, 13) ' ) ' NODE = I 

WRITE(NCOM, ' ( A, El 6 . 9 ) ' ) ' T = , TNO(I) 

45 CONTINUE 


45 



30 CONTINUE 

WRITE (NCOM, ' (A) ' ) 
GOTO 5 

100 WRITE(NCOM, ' (A) ' ) 
1200 RETURN 
END 


END DEFINE^ TEKPERATURE ' 


'♦END' 
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Appendix III Integrated Thermal / Structural Analysis 


This Appendix shows the COMET runstream of an example of integrated 
thermal/structural analysis in COMET. The SPAR element was used in thermal analysis 
and ESI element was applied in structural analysis. Different meshes were used in thermal 
and structural analyses as shown in Fig. 2. The results of a temperature distribution, a 
stress distribution and the plate buckling sensitivity parameter are also shown in Fig. 2. 
Detailed results are published in the proceedings of 35th SDM conference (AIAA Paper 
No. 94-1594). 

COMET subroutine temptran.ams which transfers temperature between different 
meshes is also listed in this appendix for future reference. 
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T 
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1 . 

Applied surface heating of magnitude q, 
BTU/in^-sec, over I by 2 region 

8 * \. 

L. 

* > 


/J\ l— 


-* 1 

x-axis is axis of symmetry 


(a) Geometry. Dimensions are in inches. 



(b) Thermal model 



Variation of thermal buckling sensitivity paramter, 
U* as a function of time for the square plate with a cutout. 



rtiom, *r 


1M«F M edftt 2-4 
W 


Temperature distribution, °F, at t cr = 13 sec for the 
square plate with a cutout. 




Cootoar Ubd 

(c) Structural model pd 

Thermal stresses a % , psi, at t Cf = 13 sec for the 
square plate with a cutout. 


Fig. 2 Thermal / structural analysis meshes and results 
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COMET Runstream 
(Thermal / Structural Analysis) 
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* rr 


cp /csm/prc/proclib.gal 
chmod +w proclib.gal 
ime testbed << \endinput 
open 28, 'proclib.gal' 

♦set plib « 28 
♦set echo*off 
♦open/new l,new.DBC 
♦add sim.clp 
♦call PLATE 
♦toe 1 

. * set echo on, ma, md 

♦open/new 2, bbb.DBC 

♦copy 2 *1,N0DA.TEMP. 1. 1 

♦copy 2 »1,DEF.K41.3 .4 

♦copy 2 , XYZ . K41 . 1 . 1 -1, JLOC.BTAB.2.5 

♦open/new 1, TEST.DBC 

♦copy 1 «2 , NODA.TEMP .1.1 

♦copy 1 *2 , DEF . K41 . 3 . 4 

♦copy 1 * 2 , XYZ . K41 . 1 . 1 

♦toe 2 

♦close 2 

♦add simply3.clp 
♦toe 1 

*def/i ES_GCP=<true> 

♦call BUCK 
♦toe 1 
[XQT EXIT 
\eof 
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♦procedure PLATE 
♦set echo-off 
♦def time-27 . 

*def/i tpl«23 
♦def/i tp2=28 
♦def temp-20. 

[XQT AUS 

♦def kx=<<226 . /3600 . >/12 . > 

TABLE (NI=9,NJ=1) :COND PROP 1:1=1 23456 
. SANDY'S Properties 

J«1 : 70 . .323 0.0906 .006122 .006122 .006122 
TABLE (NI-8 ,NJ=1) :COND COEF 1:1=1 234 
J-l: 0.323 0.0 0.0 0.0 
TABLE (NI-1, NJ-1) : KTHIC:J«1:1. 

TABLE (NI«1,NJ*1) : K AREA : J-l: 1-0 

. TABLE (NI=1,NJ=32) : SOUR K21 1 
TABLE (NI-1, NJ-16 ) : SOUR K41 2 
j-l, 16:5.0 
. J-l, 4: 0. 

. J-5 , 8 : 10.0 
. J-9,12: 0.0 
. J-13 , 16 : 10.0 
. J-17,20: 0.0 
. J-21,24: 10.0 
. J-25 , 28 : 0.0 
. J-29 , 32 : 10.0 

TABLE (NI=1,NJ=49) :TEMP NODE: DDATA-9 . : J=1 ,9:9. 

DDATA-8. :J*10, 13:125. 
DDATA-8 . : J-14 ,17:193. 
DDATA-1 . : J-18 ,24:267. 
DDATA-7. :J»25, 31:224. 
DDATA-1. :J-32, 38:210. 
J-39 : 182 . 

TABLE (NI«1,NJ=17) :TEMP NODE: DDATA-1 .: J=1 , 7 : 267 . 

DDATA-7. :J=4 0,46: 218 
DDATA-4 . : J=47 , 49 : 150 
TABLE (NI = 1,NJ=4 9) :APPL TEMP: J-l, 9 : 50 . 0 

J-10, 13:50.0 
J-14, 17: 50.0 
J-18 , 24 : 50 . 0 
J-25, 31:50.0 
J-32, 38:50.0 
J-39: 50.0 

. TABLE (NI-1 , NJ-17) :APPL TEMP : J-l , 7 : 50 . 0 

J-40, 46:50.0 
J-47, 49:50.0 

♦call MESH 
[XQT TGEO 
[XQT TRTB 

RESET Tl-0. T2=<time> DT=0.1 $ BETA-0.666667 
RESET mxndt=5000 
TEMP-0 . 0 
TSAVE-1 . 0 
[xqt aus 

define t = l tran temp 1 1 <tpl> <tp2> 
noda temp=union {t ) 

[XQT DCU 

PRINT 1 SOUR K21 
PRINT 1 NODA TEMP 1 1 
TOC 1 
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* en d 

♦procedure MESH 

.Initial Model Def. Procedure 


[xqt TAB 
START 273 
JLOC 
1 14. 0. 
9 14. 2. 
82 12. 0 
4 12. 4. 
118 14. 

8 12.5 4 
150 8. 0 
4 8. 8. 
186 12. 

8 9.8. 
218 7. : 
7 0. 2. 


0. 

16. 0. 0 

0. 

16. 2. 0 

. 0 

. 13.5 0. 

0. 

13.5 2.5 

2.5 

0. 16. 2 

. 0 

. 16. 4. 

1. 0 

. 11. o. 

0. 

11. 5. 0. 

5. 

0. 16. 5. 

0. 

16. 8. 0. 


9 19 


0 

0 

.5 

0 . 


4 19 


8 14 


4 19 


8 14 


0. 7. 8. 0 
0 . 8 . 0 . 


7 18 


[xqt ELD 

RESET NOTED =1 
K41 

MATERIAL=PROP 
NMAT*1 
NSECT - 1 
GROUP 1 

1 2 11 10 1 4 8 
GROUP 2 

5 6 15 14 1 4 4 
GROUP 3 

41 42 51 50 1 4 4 
82 83 87 86 1 3 8 
85 1 10 89 
89 10 19 93 
93 19 28 97 
97 28 37 101 
101 37 46 105 
105 46 55 109 
109 55 64 113 
113 64 73 117 
73 74 118 117 181 
117 118 126 116 

118 119 127 126 173 

116 126 134 115 

115 134 142 114 
150 151 155 154 1 3 8 
153 82 86 157 
157 86 90 161 
161 90 94 165 
165 94 98 169 
169 98 102 173 
173 102 106 177 
177 106 110 181 
181 110 114 185 
114 142 186 185 
142 143 187 186 171 

185 186 194 184 
184 194 202 183 
183 202 210 182 

186 187 195 194 1 7 3 
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158 162 219 218 
162 166 220 219 
166 170 221 220 
170 174 222 221 
174 178 223 222 
178 182 224 223 
218 219 226 225 167 

*end 
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♦procedure BUCK 

. Initial Model Def. Procedure 

♦call ES (function - 'DEFINE ELEMENTS'; esjproc * ESI; es_name » EX47) 
[xqt TAB 
START 361 6 
JLOC 


1 8. 0. 0. 

10. 0. 0. 7 1 9 

7 8. 2. 0, 

10. 4. 0. 

64 7.75 2. 

0. 9.5 4. 0. 7 1 12 

7 5. 2. 0. 

5. 4. 0. 

148 11. 0. 

0. 16. 0. 0. 619 

6 11. 5. 0 

. 16. 8. 0. 

202 3. 2. 

0. 4.5 2. 0. 4 1 7 

4 3. 4. 0. 

4.5 4. 0. 

230 2.5 5. 

0. 10.5 5. 0. 16 1 

16 0. 8. 0 

. 15. 8. 0. 

326 2.5 2. 

0. 2.5 4. 0. 6 1 6 

6 0. 2. 0. 

0. 7. 0. 


[xqt GCP 

fabrication 

shell 

fabid * 1 
nlayers = 1 
matid * 1 
thick « 0.1 
intpts * 2 
angle * 0.0 
end 
endfab 
material 

isoel 111 

1.0E7 0.3 0.322 1.0E-6 0.0 0.0 0.0 
endmat 

[xqt ELD 
<ES_EXPE CMD> 

NSECT - 1 


1 2 

9 6 

1 6 

20 

148 

149 

155 

154 

7 148 154 14 

14 

154 

160 

21 

21 

160 

166 

28 

28 

166 

172 

35 

35 

172 

178 

42 

42 

178 

184 

49 

49 

184 

190 

56 

56 

190 

196 

63 

63 

196 

245 

70 

70 

245 

244 

77 

77 

244 

243 

84 

84 

243 

242 

91 

91 

242 

241 

98 

98 

241 

240 

105 


105 

240 

239 

112 

112 

239 

238 

119 

119 

238 

237 

126 

126 

237 

236 

133 

133 

236 

235 

140 

140 

235 

234 

147 
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230 231 247 246 1 15 5 

245 196 197 261 

261 197 198 277 

277 198 199 293 

| v 293 199 200 309 

i 309 200 201 325 

202 203 207 206 136 

205 141 142 209 

F ' 209 142 143 213 

fci: 213 143 144 217 

217 144 145 221 

221 145 146 225 

| 225 146 147 229 

226 227 231 230 131 
229 147 234 .233 
326 202 206 327 
! 327 206 210 328 

1 328 210 214 329 

329 214 218 330 

330 218 222 331 

I 326 327 333 332 155 

226 230 331 222 

331 230 246 337 
337 246 262 343 
343 262 278 349 
349 278 294 355 
355 294 310 361 

[xqt TAB 
CON 1 

ZERO 6: 1,289 

i ZERO 1234: 356,361,1 

ZERO 1235: 311,325,1: 

ZERO 12345: 310: 201 
ZERO 1234: 153,195,6 
ZERO 2 4: 1,7,1 
ZERO 2 4: 148,152,1 

stop 

[XQT temptran3 
♦open 1,TEST.DBC 
♦toe 1 
[xqt ESI 

. *set echo on, ma, md 
*def/a es_name » EX47 
i *do $k » 1,6 

*def/i itime*<$k> 

*def/i ES_LOAD — SET » < i t ime > 

DEFINE TEMPERATURE 
♦find 1, DEF.EX47 . ♦ /seq=dsn 
♦rem DEF.EX47.* : dsn=<dsn> 

♦find 1, NODE. TEMP /seq=dsntemp 
♦rem NODE. TEMP : dsnjloc=<dsntemp> 

. *f ind 1, TRAN. TEMP ?seq=dsntemp 
. *rem TRAN. TEMP : dsnjloc=<dsntemp> 

♦do $j«l, 324 • 324 Elements 

ELEMENT « <$j> 

; *do $i * 12, 15 . Node #1 - #4 in DEF.EX47. 

*def/i ioffl = <$i> 

*gal2mac /Name=nnode /Iof f =<iof f 1> -- 
V'-’; /Maxn=l 1, <dsn> DATA.<$j> 
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NODE - <<$i>-ll> 

*def/i iof £2 - <<nnode>-l> J 

*gal2mac /Name-temp /Iof f *<iof f 2> /Maxn=l 1, <dsntemp> 
DATA. <<itime>> 

T=<temp> 

T=<<temp>-70 . 0> 

♦enddo 

♦enddo 

ENDJDEF INE_TEMPERATURE 
♦enddo 


[xqt DCU 

. print 1 MATL.DATA 

stop 


♦toe 1 

♦do $i»6,6 

♦ call L__ STATIC (location=centroids ; nval_metn»l 
print*<true>; load_set»<$i>) 
♦call L_STABIL_2 (function* ' eigen' ; cons_set*l; 

print=<true> ; load_set*<$i>) 

♦toe 1 
♦enddo 
♦end 


stress=<true>; 
boon set=l; n modes=3; 
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Subroutine TRANS 

(Transfer temperature between different meshes) 
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u o 


C$Header: /csm/sam/mods/RCS/temptran.ams, 93/09 maring Exp $ 
C=DECK TRANS 
O BLOCK FORTRAN 
COAUTHOR David Xue 

C=DATE September, 93 .... . 

C= PURPOSE Transfer node tempretures between different meshes 

C«IF EXTP 

PROGRAM TEMPMESH 
CALL TRANS 
STOP 
END 


c 

cc 

c 

c 

C 

C 


C=ENDIF 

c Nl-NOD^NUTOE^OF^SHl, N2=NODE NUMBER OF MESH2 , N3 .ELEMENT NUMBER OF MESH2 
c Transfer temperature from mesh2 to meshl 

PARAMETER (Nl*361 , N2=273 , N3=240 , NELN=4 ) 

COMMON/MAT/A1 (3,N1) , A2 (3,N2) ,A3 (N3, 8) ,NA4 (4,N3) , INODEN(Nl) , 

$ TA1 (Nl, 6) , TA2 (N2 , 6) 

DIMENSION NA5(16,N3) 

CHARACTER* 4 0 dsl , ds2 , ds3 , ds4 , rnl , m2 , m3 , rn4 , rn5 
Read temperature from NODA . TEMP . 1 . 1 (6 records) to TA2(N2,6) 
lib=l 

lib=LMOPEN {'COLD', 0, ' TEST . DBC' , 0, 2000) 

CALL GMCODN (dsl, 'NODA ','TEMP ',1,1,0) 
iesn=LMFIND (lib, dsl, 1800) 

CALL GMCORN (rnl, 'DATA ',1,6) 

CALL GMGETN ( ' R' , lib, iesn, rnl , ' S' , TA2 ,N, 0,0, 0,1800) 

CALL GMGETN ('R' , lib, iesn, rnl , 'S' ,TA2, -N2 ,N2 , 0, 0, 1801) 
print *, 'TA2(N2,6)«' 
do 30 1*1, N2 

30 print *,'N2*',I,' * ' , (TA2 (I , J) , J*1 , 6 ) . . , - 

Read node coordinates o£ K41 from XYZ.K41.1.1 which is copied from old 
JLOC . BTAB .2.5 of K41 element meshes. 

CALL GMCODN (ds2 , ' XYZ ','K41 ',1,1,0) 
iesn2*LMFIND (lib, ds2 ,1700) 

CALL GMCORN (m2, 'DATA ',1,0) 

CALL GMGETN ( ' R' , lib, iesn2 , m2 , 'S' , A2,N, 0, 0, 0, 1701) 
print *, ' OLD K41 MESH NODE COORD. X, Y, Z, A2(3,N2)*' 
do 35 1=1, N2 

35 print *,'N2=',I,' * ' , (A2 ( J, I) , J=1 , 3) 

C Read node coordinates of EX47 from JLOC . BTAB . 2 . 5 
CALL GMCODN (ds 3, 'JLOC ' , ' BTAB ',2,5,0) 
iesn3=LMFIND (lib, ds3 ,1600) 

CALL GMCORN (m3 , ' DATA ',1,0) 

CALL GMGETN ('R' ,lib, iesn3,m3, 'S' ,A1,N, 0, 0, 0, 1601) 
print * , 'NEW EX47 MESH NODE COORD. X, Y, Z, A1(3,N1)=' 
do 33 1*1, Nl 

33 print *,'N1*',I,' *' , (A1 (J, I) , J*l,3) 

C Read element-node numbers of K41 from DEF.K41.3.4 
CALL GMCODN (ds4, 'DEF ','K41 ',3,4,0) 
iesn4*LMFIND (lib, ds4, 1500) 

CALL GMCORN (rn4, 'DATA ',1,N3) 

CALL GMGETN ( 'R' , lib, iesn4 , rn4 , ' I' ,NA5,N, 0,0,0,1501) 

DO 41 J1=1,N3 
DO 41 J2 = 1 , 4 

41 NA4 ( J2 , Jl) *NA5 ( J2 + 12 , Jl) 

c CALL GMGETN ('R/L' , lib, iesn4 , rn4 , ' I' ,NA4 , 16 , 0 , 0, 12, 1501) 


print * , ' NA5 . 1* ' , (NA5 ( J, 1) , J=1 , 16 ) 
print * , ' NA5 . 2= ' , (NA5 ( J, 2 ) , J=1 , 16 ) 


58 



o o 


print * , ' NA4 (4 , N3 ) - ELEMENT_NODES ' 
do 45 I-1,N3 

45 print * , ' ELEMENT* ' , I , ' * ' , (NA4 ( J, I) , J*1 , 4 ) 

CALL TTT 

C Create new temperature NODE. TEMP for EX47 meshes 
CALL GMPUNT( lib, 'NODE. TEMP idsn, 4800, 2500) 

c CALL GMCORN(rn5, 'DATA ',1,1) 

CALL GMCORN(rn5, 'DATA ',1,6) 

c CALL GMPUTN (' W' , 1, idsn, rnS, ' S' ,TA1, -Nl, 1,0, 0,2501) 

CALL GMPUTN ('W' , 1, idsn, rn5, ' S' ,TA1, -Nl, 6, 0, 0, 2501) 

CALL GMCLOS (lib, 0,2600) 

c lib=LMOPEN ('COLD', 0, ' TEST. DBC' , 0, 2000) 

RETURN 

END 

C=END FORTRAN 
c 

C-DECK TTT 
C=BLOCK FORTRAN 

SUBROUTINE TTT 

c PARAMETER (Nl*121 , N2*81 , N3=64 , NELN=4 ) 

PARAMETER (Nl=3 6 1 , N2 «2 73 , N3 *24 0 , NELN=4 ) 

COMMON /MAT/ A1 (3,N1) ,A2 (3,N2) ,A3 (N3, 8) ,NA4 (4,N3) , INODEN (Nl) , 
$ TA1 (Nl, 6) , TA2 (N2, 6) 

DIMENSION EXMAX (N3) ,EXMIN(N3) , EYMAX (N3) , EYMIN (N3 ) 

C A3 (N3 , 8 ) =xl , x2 , x3 , x4 , yl , y2 , y3 , y4 
DO 2 1*1, N3 
DO 5 J* 1 , NELN 
A3 ( I , J) *A2 ( 1 , NA4 ( J, I ) ) 

A3 (I, J+NELN) *A2 (2,NA4 (J,I) ) 

5 CONTINUE 

2 CONTINUE 

print *, ' A3 (N3, 8) * ELEMENT COORDINATES' 
do 7 1*1, N3 
c 7 print *,'I«',I,' *' , (A3 (I, J) , J*l, 8) 

DO 10 1*1 , Nl 
10 INODEN ( I) *0 

NTOTAL-O 

C Find all common nodes both meshes have 
DO 20 1*1, Nl 
X1=A1 (1,1) 

Y1=A1 (2, 1) 

Z1*A1 (3,1) 

DO 25 J-1,N2 
X2=A 2 (1, J) 

Y2*A 2 (2, J) 

Z2«A2<3,J) 

XR1*ABS (X1-X2) 

XR2*0 . 01*X1 

IF (XR1 . GT . XR2) GOTO 25 

YR1*ABS (Y1-Y2) 

YR2=0 . 01*Y1 

I F ( YR1 . GT . YR2 ) GOTO 25 

DO 27 K=1 , 6 

TA1 ( I , K) =TA2 ( J, K) 

c print *, 'K*' ,K, ' *1*' , I, ' *J=' , J, ' TA1 { I , K) * ' , TA1 ( I , K) 

27 CONTINUE 

C INODEN (I) *1, Ith node temperature has been evalueted. 

INODEN (I) =1 
NTOTAL=NTOTAL+ 1 
IF (NTOTAL.EQ.N1) GOTO 100 
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25 CONTINUE 

20 CONTINUE 

print *, ' Numbers of 
C Find Xmax, Xmin, Ymax and 
DO 30 1*1, N3 
XMAX*0 .0 
XMIN-1.0E20 


equare nodes are: ',NTOTAL 
Ymin for each element 


XMIN=X 

YMAX«Y 

YMIN*Y 


XMAXMIN, YMAXMIN= ' , XMAX, XMIN, YMAX, YMIN 


YMAX*0 . 0 
YMIN*1 . 0E20 
DO 35 J*1,NELN 
X=A3 (I, J) 

Y=A3 ( I , NELN+ J) 

IF (XMAX.LT.X) XMAX»X 
IF (XMIN.GT.X) 

IF (YMAX.LT. Y) 

IF {YMIN.GT. Y} 

35 CONTINUE 

EXMAX ( I ) =XMAX 
EXMIN (I) =XMIN 
E YMAX (I) =YMAX 
EYMIN(I) “YMIN 

c print ELEMENT* ', I 

30 CONTINUE 

C Find if Jth node is in Ith element 
DO 40 1*1, N3 

AX1»A3 (1,1) 

AX2-A3 (1,2) 

AX3=A3 (1,3) 

AX4 *A3 (1,4) 

AY1*A3 (1,5) 

AY2*A3 (1,6) 

AY3-A3 (1,7) 

AY4 »A3 (1,8) 

CALL AAR (AREA12 , AX1 , AX2 , AX3 , AY1 , AY2 , AY3 ) 

CALL AAR ( AREA3 4 , AX1 , AX3 , AX4 , AY1 , AY3 , AY4 ) 

AREA® AREA1 2 +AREA3 4 
XMAX “EXMAX (I) 

XMIN»EXMIN ( I ) 

YMAX =E YMAX ( I ) 

YMIN“EYMIN ( I ) 

DO 45 J»1,N1 

IF (INODEN ( J) . EQ . 1) GOTO 45 
X“A1(1, J) 

IF(X GT. XMAX. OR. X.LT. XMIN. OR. Y.GT. YMAX. OR. Y.LT. YMIN) GOTO 45 

CALL AAR ( AREA1 , X , AX1 , AX2 , Y , AY1 , AY2 ) 

IF (AREA1.GT. AREA) GOTO 45 

CALL AAR ( AREA2 , X, AX2 , AX3 , Y , AY2 , AY3 ) 

AX 1 2 “AREA1 +AREA2 

IF (AX12 .GT . AREA) GOTO 45 

r&T.T. RAR ( AREA3 . X . AX3 . AX4 , Y , AY3 , AY4 ) 


AX12 3=AX12 +AREA3 
IF (AX123 . GT . AREA) GOTO 45 
CALL AAR (AREA4 , X, AX4 , AX1 , Y , AY4 , AY1) 
AX1 2 3 4 =AX1 2 3 +AREA4 
ERROR=ABS { (AX1234 -AREA) /AREA) 
c PRINT * , ' ERROR* ' , ERROR 

IF (ERROR. GT. 0 . 01) GOTO 45 

C Find the Jth node is in triangle 1-2-3 or 3-4-1 
IF (AX12 .GT. AREA12) THEN 
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CALL 


50 


53 

54 


45 

40 

100 


DO 50 K-1,6 
T1-TA2 (NA4(1,I),K) 

T2-TA2 (NA4 (2,1) , K) 

T3-TA2 (NA4 (3, I) , K) 

T4-TA2 (NA4 (4, I) , K) 

CALL INTEP (AREA34 , T3 , T4 , T1 , AX3 , AX4 , AX1 , AY 3 , AY4 , AY1 , TEMP , X , Y) 
TA1(J,K)*TEMP 
CONTINUE 
GOTO 54 
END IF 

DO 53 K*l,6 
T1-TA2 (NA4 (1, I) ,K) 

T2*TA2 (NA4 (2, 1) , K) 

T3=TA2 (NA4 (3,1) ,K) 

CALL INTEP^AREA12?T1 ! «! T3 , AX1 , AX2 , AX3 , AY1 , AY2 , AY3 , TEMP, X , Y) 
TA1(J,K)*TEMP 
CONTINUE 

INODEN(J) *1 

NTOTAL=NTOTAL+ 1 , J . _ T 

print * , ' new node ',J,' is m old element ,1 
IF(NTOTAL.EQ.Nl)GOTO 100 
CONTINUE 
CONTINUE 

CONTINUE 


c **************************************** 
c FOR PATRAN INPUT RESULT FILE (SEE PATRAN MANUAL 27.7.2.1) 

NNODES1-N1 

NNODES2=N2 

MAXNOD-O 

DEFMAX=0 . 0 

NDMAX*0 

NWIDTH-1 

Q **************************************** 

DO 70 K=6 , 6 

* t ' K * ' ,K, ' * OLD TEMPERATURE IS 
print * , (TA2 ( I , K) , 1*1 , N2 ) 
print *,'K«',K, ' * NEW TEMPERATURE IS' 
print * , (TA1 ( I , K) , 1*1 , Nl) 

70 CONTINUE , 

open (unit=9, file- 'TEMP. NEW' , status=' old' ) 
write (9, *) 'TEMPERATURE PATRAN RESULT FILE' 

write (9, 7111) NNODES1, MAXNOD, DEFMAX, NDMAX, NWIDTH 

7111 FORMAT (219, E15. 6, 219) 

write ( 9 , * ) ' IN PATRAN FORMAT ' 
write (9, *) 'NEW NODE TEMPERATURE' 

DO 72 1*1, N1 

write (9,7222) I , TA1 (1,6) 

7222 FORMAT ( I 8 , El 3 . 7 ) 

72 CONTINUE 

close (unit*9) 

open (unit=9, file* 'TEMP .OLD' , status*' old' ) 
write (9, *) 'TEMPERATURE PATRAN RESULT FILE' 

write (9, 7 111 )NNODES 2 , MAXNOD, DEFMAX, NDMAX, NWIDTH 
write (9, *) ' IN PATRAN FORMAT' 
write (9, *) 'OLD NODE TEMPERATURE' 

DO 74 1=1, N2 

write (9,7222) I , TA2 (1,6) 

74 CONTINUE 

close (unit=9) 
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C=END FORTRAN 


return 

end 


ODECK AAR 
C=BLOCK FORTRAN 

SUBROUTINE AAR(A,X1,X2,X3,Y1, Y2,Y3) 

A*0. 5* ( (X3-X1) * (Y1-Y2) - (X1-X2) * (Y3-Y1) ) 
A=ABS (A) 

RETURN 

END 

C=END FORTRAN 


C=DECK INTEP 

C=BLOCK FORTRAN _ . 

SUBROUTINE INTEP (A, T1 , T2 , T3 , XI , X2 , X3 , Y1 , Y2 , Y3 , . 

A1«Y2-Y3 

B1»X3-X2 

C1=X2*Y3 -Y2*X3 

XN1- (A1*X+B1*Y+C1) *0 . 5/A 

A2=Y3 -Y1 


B2=X1 -X3 

C2*=X3 *Y1 -Y3 *X1 

XN2= (A2*X+B2*Y+C2) *0 . 5/A 

A3*Y1-Y2 


B3=X2 -XI 
C3=X1*Y2-Y1*X2 
XN3- (A3*X+B3*Y+C3) *0 . 5/A 
T«XN1*T1+XN2*T2+XN3*T3 
RETURN 
END 

C-END FORTRAN 


,X,Y) 


62 



